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O ! Abstract 
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^ , It is argued that, for motion in a central force field, polar reciprocals of trajectories are an elegant 

alternative to hodographs. The principal advantage of polar reciprocals is that the transformation 
from a trajectory to its polar reciprocal is its own inverse. The form of polar reciprocals A;* of 
Q^, Kepler problem orbits is established, and then the orbits k themselves are shown to be conic 

lyj ' sections using the fact that k is the polar reciprocal of fc*. A geometrical construction is presented 

'T^ ■ for the orbits of the Kepler problem starting from their polar reciprocals. No obscure knowledge 

Y^ . of conies is required to demonstrate the validity of the method. Unlike a graphical procedure 

£^ ' suggested by Feynman (and amended by Derbes), the algorithm based on polar reciprocals works 
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without alteration for all three kinds of trajectories in the Kepler problem (elliptical, parabolic, 
and hyperbolic). 



Approximately, a hodograph is a plot of velocities along a trajectory; more precisely, it 
is the locus of the tips of the velocity vectors after they have been parallelly transported 
until their tails are at the origin (in velocity space). There have been many articles on the 
pedagogic virtues of hodographs.^"^ Feynman's Lost Lecture'^ contains, amongst other things, 
a recipe for drawing the elliptical path of a planet given its hodograph. The procedure, as 
reproduced in Ref. 6, has some shortcomings,^'^ but, fortunately, these have been more 
than satisfactorily rectified by Derbes.^ Derbes notes that Feynman's scheme also works for 
hyperbolic orbits and he is able to devise another construction for the exceptional case of 
parabolic orbits. A completely different way of tracing all three kinds of orbits with the help 
of hodographs has been successfully developed by Salas-Brito and co-workers in a series of 
publications culminating in Ref. 10. 

The various geometrical methods of the previous paragraph are easily implemented, but, 
to prove their validity, a student would have to be acquainted with many properties of conies 
which no longer form part of most school curricula. Indeed, at least one reader of Derbes' 
article feels that his exposure to the hodograph left him "disappointed by the trade-off of 
intricate calculus for obscure geometry". ^^ There is, however, an alternative approach which 
requires only some elementary vector algebra and calculus for its justification. 

The seed for this other construction can be traced back to a result of Newton (Proposition 
I, Corollary I on page 41 of Ref. 12) arising from the conservation of angular momentum 
in a central force field. Let r be the position vector of a body relative to the center O of 
the force field and let C be the body's angular momentum per unit mass with respect to O; 
then, 

C = |C| = |r X r I = r_L|r I , (1) 

where r_|_ is the component of r perpendicular to the instantaneous velocity v = r (relative 
to O): since C is a constant of the motion, the speed v = |r | of the body is inversely 
proportional to r_|_ (z/ = C/r^). Consider now the mapping P ^ P^ illustrated in Fig. 1. 
The point P is on a typical trajectory r of a body in an (attractive) central force field (with 
center O); of course, r is in the plane through O perpendicular to C (and so are r and u). 
By construction, the vector OP^, (in this plane) is perpendicular to the velocity z/ at P and 
has magnitude |0P*| = r^^. Since rj^ = u/C [from Eq. (1)] and the direction of OP^ is 
obtained from that of z/ by a clockwise rotation through 90° (see Fig. 1), the curve r* in 
Fig. 1, which is the locus of points P^ as P varies over r, is the hodograph of t, rescaled by 




FIG. 1. The map P ^>- P>f and the polar reciprocal r* of r; r* = r_|_ . (The notation r* for OP* is 
convenient in Appendix A.) 

a factor of C^^ and rotated clockwise through 9(F (see also section 3 of Ref. 13). Following 
Ref. 13, I will term r* the polar reciprocal of r. 

It is the inverse of the mapping P — ?> P* in Fig. 1 which can be used to draw trajectories. 
In fact, there is a pleasing symmetry. The mapping P — ?> P* is involutory, i.e. (P*)* = P 
(see Appendix A for an elementary proof), so trajectories in a central force field and their 
hodographs (after rotation and rescaling as in the preceding paragraph) are polar reciprocals 
of each other. The point-by-point determination of a trajectory from a rotated and rescaled 
hodograph involves exactly the steps depicted in Fig. l}^ Unlike the methods of Refs. 9 and 
10, polar reciprocation is valid for any smooth hodograph associated with any central force 
field. Be that as it may, I will now specialize to orbits of the Kepler problem. I will also use 
the identity 

OP, = vx C/C^, (2) 

_i 
which takes advantage of the fact that, for the trajectory r depicted in Fig. 1, C points 

perpendicularly out of the page, so that z/ x C is parallel to OP^ and has magnitude vC. 

For an inverse-square force per unit mass of —'-fr/r^ (7 > 0), the equation of motion for 

the position vector OP^ of a typical point P* on the polar reciprocal of an orbit (see Fig. 1) 

reads 

d -— ^ dv C 7 ^ C '-) d(j) ^ dcr , . 



where p = 7/C^ and plane polar coordinates r and have been adopted for the orbital 
plane; the corresponding unit vectors are e^ and e^ with the origin of the coordinate system 
at the force center O, and the azimuthal angle defined as in Fig. 2 (so that C = r'^(j)k). 
Equation (3) implies that OP^, — pe,,. is a (vectorial) constant of the motion, say 00 ^^ (drawn 
in Fig. 2). Setting 00 ^^ = pe*, 

OP, = p(e, + er), (4) 

which demonstrates that the polar reciprocal k^ of an orbit k of the Kepler problem is a 
circle or an arc of a circle (as suggested by the plot in Fig. 2).^^ 

The polar reciprocal of A;* (which would be the corresponding orbit k) comprises points 
like P** in Fig. 2, the image (under polar reciprocation) of P* on /c*. In terms of the angles 
in Fig. 2, r = OP^^: = l/OQ = [OP* cos(0 — 0*)]~^ = {OP^ ■ e^)"^ or, evaluating OP^: ■ e^ 
with Eq. (4), 

- = p(e* ■ Cr + 1) = p(l + |e*| COS0). (5) 

Comparison of Eq. (5) with the standard equation^^ for a conic in polar coordinates confirms 
that P** is on a conic with focus O, eccentricity e = |e*|, latus rectum 2/p, and directrix 
perpendicular to 00 ^:. The angle can thus be identified as the true anomaly (of celestial 
mechanics), and e* can be reinterpreted as a vector of magnitude equal to the eccentricity 
e of the orbit k, directed from the force center O to the point on k of closest approach (i.e. 
the periapse). In fact, to within a constant, e* is the Laplace- Runge-Lenz vector. ^^ 

In the previous two paragraphs, it has been shown that it is easy to establish the character 



FIG. 2. A typical point P** on the polar reciprocal of A;*; ZP^,OP^:.^ = — 0* as O-^P^: \\ OP^.^. By 
construction, OP** = OQ-^ = [OP^cosiPX)K*)r^ just as, in Fig. 1, OP^ = [OPcos(POK)]"^. 



of the polar reciprocal A;* of an orbit k of the Kepler problem and even easier to infer from /c* 
that the orbit k must be a conic section. It is now possible to indicate how, as an alternative 
to the methods of Refs. 9 and 10, polar reciprocation may be used, in principle, to draw 
orbits with a compass and ruler. 

Suppose that the orbiting body's velocity is given at some point Pq on the orbit (which 
need not coincide with an apse of the orbit); let tq be the position vector of Pq relative to 
O, and let the associated velocity (in the rest frame of O) be uq. First, one must construct 
/c*. The following sequence of steps can be employed (Pq* below denotes the image under 
polar reciprocation of Pq): 

(i) calculate C = {tq x i)q\ and then p (= '~fC~'^); 



(ii) add graphically the vectors OPq* = ^'o x C/C'^ and Pq^O^ = — (p/ro)ro to find the 
position of O* relative to O; 

(iii) determine the eccentricity e (= 00^/ p); 

(iv) at O*, draw a circle of radius p if e < 1 or a circular arc on which the points P* are 
such that |ZOO*P*| > cos^^ e^^ if e > 1 (cf. Fig. 3 in Appendix B).^^ 

From a point P* on /c*, the corresponding point P** on k can be determined by demanding 
that OP^:^: is parallel to 0^:P^ and has magnitude |0P**| = 1/OQ. Well-known compass- 
and-ruler constructions suffice. One can: 

(1) erect at P* a fine (/_l) perpendicular to the line through O^^P^,-^ 

(2) draw through O the line (/j|) parallel to the line segment O^^P^,, and then; 

(3) invert the point of intersection of l± and /y (the point Q in Fig. 2) with respect to the 
unit circle centered on O. 

Compass-and-ruler implementations (with JavaSketchpad) of all three of these procedures 
can be found at <http://www.susqu.edu/brakke/constructions/constructions.htm> 
(accessed July 2, 2011).^^ 

For students today, who have access to programs like Mathematica, it may seem that 
there is little or no need for the geometrical construction under discussion; it takes only a 
few keystrokes (as in Example 2.1.9 of Ref. 20) to generate a plot of a Kepler problem orbit; I 



have, nevertheless, witnessed the intellectual satisfaction students experience in being able to 
correctly anticipate the character and orientation of an orbit: these features can be deduced 
from initial conditions via steps (i) to (iii) of the preceding paragraph. A more traditional 
alternative would be to use the Laplace-Runge-Lenz vector, but I have found that even my 
most talented and dedicated junior-level students are somewhat mystified by this construct 
and disinclined to use it, when it is introduced in isolation (as it usually is) as a single- 
valued combination of dynamical variables which happens to be a constant of the motion. 
My students have been more receptive to this "exotic" constant of the motion when it is 
presented within the context of a discussion of hodographs or polar reciprocals of the Kepler 
problem. Other more mundane issues can be tackled. Some suggestions for problems are 
given in Appendix B. 
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Appendix A: PROOF THAT P** = P IN A CENTRAL FORCE FIELD 

The proof below involves some elementary vector algebra and, crucially, use of the equa- 
tions of motion, which can be assumed to be of the form 

r = f{r)r, (Al) 

where r(t) = x(t)T+ y{t)'j~ is the position of the orbiting body relative to the force center 

(the orbital plane is taken the xy-plane with the 2;-axis parallel to the angular momentum 

per unit mass C). 

By definition, under polar reciprocation, the point P with position vector r(t) (on trajec- 

tory r) is mapped to the point P* (on r^,) with position vector r*(t) = r(t) x C/C^, where 

_\ ^ ^ 

C = r{t) X r(t). Likewise, under polar reciprocation, P^, is mapped to the point P** with 

position vector r^,* = r*(t) x C^/Cl, where C* = r*(t) x r*(t). 

Substituting r(t) = x(t)T+ y{t)Jmio the definition of r^, one finds that 

r,(t) = [y(t)?-x(t)J]/C, (A2) 



where C = x{t)y{t) — y{t)x{t) (> 0). Using Eq. (Al) and the fact that C is a constant of 
the motion, r* can be simphfied to 

7^,{t) = f{r)[y{t)t-x{t)j]/C, (A3) 

and, hence, 

C, = -f{r)k/C. (A4) 

Equations (A3) and (A4) imply that r^{t) x C* = C^r or r*^, = r (as required). 

Appendix B: SUGGESTED PROBLEMS 

(1) It is not unusual to see the orbital plane (or xy-plane) identified with the complex 
plane (or 2;-plane) via the (obvious) isomorphism {x,y) -^ z = x + iy. Formulate polar 
reciprocation as an operation in the complex plane. Use this representation to prove polar 
reciprocation is involutory.^^ 

(2) The polar reciprocal k^ of a known orbit k of the Kepler problem can be easily found. 

a) Justify the claim that the images under polar reciprocation of just two points on k suffice 
to fix A;*. 

b) Show that, with the a;-axis aligned along e^,, the velocity at points around k is 

u=ilp)He?+e^) (Bl) 

in the notation of this paper. 

c) By considering the periapse A oi k (where va J- rA and rA \\ ^*), demonstrate that the 
energy E of the orbiting body is related to the eccentricity e of A; by 

- = \ip{e' - 1), (B2) 

where m is the mass of the orbiting body. (Begin by finding expressions for r^ = va/C 
and va in terms of p and e.) 

d) Use Eqs. (Bl) and (B2) to prove that the inequality z/^ (= 2E /m + 27/r) > 2E/m 
reduces to ecos0 > — 1. 

(3) The polar reciprocal of an elliptical orbit (of the Kepler problem) must be a complete 




FIG. 3. The polar reciprocal fc* of the hyperbolic orbit k of Fig. 4. The point 5* is the image 
under polar reciprocation of point B in Fig. 4; P^°° are the images of the two points at infinity 
P = lim^^-i-oo Pit) on the in/out asymptotes. Like k, k^, is symmetric about the principal axis 
of k (the line through O and O*). 

circle (as opposed to only an arc of a circle) because of the periodicity of the motion. Show 
that this circle has radius 

where r^ and r^ are the distances (from the force center O) of the periapse and apoapse, 
respectively, and that the distance of its center O* from O is 

Infer from Eqs. (B3) and (B4) an expression for the eccentricity of the orbit. 

(4) Figure 3 depicts the polar reciprocal /c* of the hyperbolic orbit k in Fig. 4. Prove that 
0P^°° is tangent to the circle in Fig. 3 of which /c^, is a part by showing that AOP^'^Otf = 
90°. By appealing to the "tangent-secant theorem", ^^ deduce that OA* = rs/h"^- Hence, 
conclude that, in Fig. 3, the circle has radius 

and 

1 / 



o''- = 2;;l^ + #'- <««) 
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FIG. 4. A hyperbolic orbit k of the Kepler problem. 

Use Eqs. (B5) and (B6) to justify the assertion that results appropriate to paraboUc orbits 
can be obtained by taking the hmit 6 — ?> oo, keeping tb fixed. 

(5) Many spacecraft maneuvers involve short firings of high-thrust rockets, which can be 
modeled as resulting in instantaneous changes Az/ in velocity with no change in position. 
The important Hohmann transfer (depicted in Fig. 5) requires two such impulsive thrusts: 
the first (at Pi) places the vehicle on the semi-elliptical trajectory hif, and the second (at 
Pf) puts the vehicle into the circular orbit cj.^^ 




FIG. 5. The Hohmann transfer trajectory hif between coplanar circular orbits Cj and c/ in 
field of a spherical mass centered on O; hif is tangent to Cj at Pi and tangent to Cf at Pf. 



the 



a) Show that, provided Az/ is tangential to the spacecraft's orbit and the direction of motion 
of the spacecraft is not reversed, r* = z/ x C jC'^ is unchanged by an impulsive thrust. 

b) Find the polar reciprocal representation of the Hohmann transfer in Fig. 5, i.e. draw (in 
one diagram) the polar reciprocals of q, C/, and /ij/. 

c) Generalize the result of part b) to the case of two counter-clockwise non-intersecting co- 
planar and co-axial elliptical orbits about the force center O. The transfer trajectory must 
connect the periapse of the inner orbit to the apoapse of the outer orbit. 
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